RNAseq based variant dataset in a black poplar association panel.

OBJECTIVE: Black poplar (Populus nigra L.) is a species native to Eurasia with a wide distribution area. It is an ecologically important species from riparian ecosystems, that is used as a parent of interspecific (P. deltoides x P. nigra) cultivated poplar hybrids. Variant detection from transcriptomics sequences of 241 P. nigra individuals, sampled in natural populations from 11 river catchments (in four European countries) is described here. These data provide new valuable resources for population structure analysis, population genomics and genome-wide association studies. DATA DESCRIPTION: We generated transcriptomics data from a mixture of young differentiating xylem and cambium tissues of 480 Populus nigra trees sampled in a common garden experiment located at Orléans (France), corresponding to 241 genotypes (2 clonal replicates per genotype, at maximum) by using RNAseq technology. We launched on the resulting sequences an in-silico pipeline that allowed us to obtain 878,957 biallelic polymorphisms without missing data. More than 99% of these positions are annotated and 98.8% are located on the 19 chromosomes of the P. trichocarpa reference genome. The raw RNAseq sequences are available at the NCBI Sequence Read Archive SPR188754 and the variant dataset at the Recherche Data Gouv repository under https://doi.org/10.15454/8DQXK5 .

p-Coumaroylation of poplar lignins impacts lignin structure and improves wood saccharification.

The enzymatic hydrolysis of cellulose into glucose, referred to as saccharification, is severely hampered by lignins. Here, we analyzed transgenic poplars (Populus tremula × Populus alba) expressing the Brachypodium (Brachypodium distachyon) p-coumaroyl-Coenzyme A monolignol transferase 1 (BdPMT1) gene driven by the Arabidopsis (Arabidopsis thaliana) Cinnamate 4-Hydroxylase (AtC4H) promoter in the wild-type (WT) line and in a line overexpressing the Arabidopsis Ferulate 5-Hydroxylase (AtF5H). BdPMT1 encodes a transferase which catalyzes the acylation of monolignols by p-coumaric acid (pCA). Several BdPMT1-OE/WT and BdPMT1-OE/AtF5H-OE lines were grown in the greenhouse, and BdPMT1 expression in xylem was confirmed by RT-PCR. Analyses of poplar stem cell walls (CWs) and of the corresponding purified dioxan lignins (DLs) revealed that BdPMT1-OE lignins were as p-coumaroylated as lignins from C3 grass straws. For some transformants, pCA levels reached 11 mg·g-1 CW and 66 mg·g-1 DL, exceeding levels in Brachypodium or wheat (Triticum aestivum) samples. This unprecedentedly high lignin p-coumaroylation affected neither poplar growth nor stem lignin content. Interestingly, p-coumaroylation of poplar lignins was not favored in BdPMT1-OE/AtF5H-OE transgenic lines despite their high frequency of syringyl units. However, lignins of all BdPMT1-OE lines were structurally modified, with an increase of terminal unit with free phenolic groups. Relative to controls, this increase argues for a reduced polymerization degree of BdPMT1-OE lignins and makes them more soluble in cold NaOH solution. The p-coumaroylation of poplar samples improved the saccharification yield of alkali-pretreated CW, demonstrating that the genetically driven p-coumaroylation of lignins is a promising strategy to make wood lignins more susceptible to alkaline treatments used during the industrial processing of lignocellulosics.

RNAi suppression of DNA methylation affects the drought stress response and genome integrity in transgenic poplar.

Trees are long-lived organisms that continuously adapt to their environments, a process in which epigenetic mechanisms are likely to play a key role. Via downregulation of the chromatin remodeler DECREASED IN DNA METHYLATION 1 (DDM1) in poplar (Populus tremula × Populus alba) RNAi lines, we examined how DNA methylation coordinates genomic and physiological responses to moderate water deficit. We compared the growth and drought response of two RNAi-ddm1 lines to wild-type (WT) trees under well-watered and water deficit/rewatering conditions, and analyzed their methylomes, transcriptomes, mobilomes and phytohormone contents in the shoot apical meristem. The RNAi-ddm1 lines were more tolerant to drought-induced cavitation but did not differ in height or stem diameter growth. About 5000 differentially methylated regions were consistently detected in both RNAi-ddm1 lines, colocalizing with 910 genes and 89 active transposable elements. Under water deficit conditions, 136 differentially expressed genes were found, including many involved in phytohormone pathways; changes in phytohormone concentrations were also detected. Finally, the combination of hypomethylation and drought led to the mobility of two transposable elements. Our findings suggest major roles for DNA methylation in regulation of genes involved in hormone-related stress responses, and the maintenance of genome integrity through repression of transposable elements.

Gene expression predictions and networks in natural populations supports the omnigenic theory.

BACKGROUND: Recent literature on the differential role of genes within networks distinguishes core from peripheral genes. If previous works have shown contrasting features between them, whether such categorization matters for phenotype prediction remains to be studied. RESULTS: We measured 17 phenotypic traits for 241 cloned genotypes from a Populus nigra collection, covering growth, phenology, chemical and physical properties. We also sequenced RNA for each genotype and built co-expression networks to define core and peripheral genes. We found that cores were more differentiated between populations than peripherals while being less variable, suggesting that they have been constrained through potentially divergent selection. We also showed that while cores were overrepresented in a subset of genes statistically selected for their capacity to predict the phenotypes (by Boruta algorithm), they did not systematically predict better than peripherals or even random genes. CONCLUSION: Our work is the first attempt to assess the importance of co-expression network connectivity in phenotype prediction. While highly connected core genes appear to be important, they do not bear enough information to systematically predict better quantitative traits than other gene sets.

Cytological, Biochemical and Molecular Events of the Embryogenic State in Douglas-fir (Pseudotsuga menziesii [Mirb.]).

Somatic embryogenesis techniques have been developed for most coniferous species, but only using very juvenile material. To extend the techniques' scope, better integrated understanding of the key biological, physiological and molecular characteristics of embryogenic state is required. Therefore, embryonal masses (EMs) and non-embryogenic calli (NECs) have been compared during proliferation at multiple levels. EMs and NECs originating from a single somatic embryo (isogenic lines) of each of three unrelated genotypes were used in the analyses, which included comparison of the lines' anatomy by transmission light microscopy, transcriptomes by RNAseq Illumina sequencing, proteomes by free-gel analysis, contents of endogenous phytohormones (indole-3-acetic acid, cytokinins and ABA) by LC-MS analysis, and soluble sugar contents by HPLC. EMs were characterized by upregulation (relative to levels in NECs) of transcripts, proteins, transcription factors and active cytokinins associated with cell differentiation accompanied by histological, carbohydrate content and genetic markers of cell division. In contrast, NECs were characterized by upregulation (relative to levels in EMs) of transcripts, proteins and products associated with responses to stimuli (ABA, degradation forms of cytokinins, phenols), oxidative stress (reactive oxygen species) and carbohydrate storage (starch). Sub-Network Enrichment Analyses that highlighted functions and interactions of transcripts and proteins that significantly differed between EMs and NECs corroborated these findings. The study shows the utility of a novel approach involving integrated multi-scale transcriptomic, proteomic, biochemical, histological and anatomical analyses to obtain insights into molecular events associated with embryogenesis and more specifically to the embryogenic state of cell in Douglas-fir.

Accuracy of RNAseq based SNP discovery and genotyping in Populusnigra.

BACKGROUD: Populus nigra is a major tree species of ecological and economic importance for which several initiatives have been set up to create genomic resources. In order to access the large number of Single Nucleotide Polymorphisms (SNPs) typically needed to carry out a genome scan, the present study aimed at evaluating RNA sequencing as a tool to discover and type SNPs in genes within natural populations of P. nigra. RESULTS: We have devised a bioinformatics pipeline to call and type SNPs from RNAseq reads and applied it to P. nigra transcriptomic data. The accuracy of the resulting RNAseq-based SNP calling and typing has been evaluated by (i) comparing their position and alleles to those previously reported in candidate genes, (ii) assessing their genotyping accuracy with respect to a previously available SNP chip and (iii) evaluating their inter-annual repeatability. We found that a combination of several callers yields a good compromise between the number of variants type and the accuracy of genotyping. We further used the resulting genotypic data to carry out basic genetic analyses whose results confirm the quality of the RNAseq-based SNP dataset. CONCLUSIONS: We demonstrated the potential and accuracy of RNAseq as an efficient way to genotype SNPs in P. nigra. These results open prospects towards the use of this technology for quantitative and population genomics studies.

Changes in the epigenome and transcriptome of the poplar shoot apical meristem in response to water availability affect preferentially hormone pathways.

The adaptive capacity of long-lived organisms such as trees to the predicted climate changes, including severe and successive drought episodes, will depend on the presence of genetic diversity and phenotypic plasticity. Here, the involvement of epigenetic mechanisms in phenotypic plasticity toward soil water availability was examined in Populus×euramericana. This work aimed at characterizing (i) the transcriptome plasticity, (ii) the genome-wide plasticity of DNA methylation, and (iii) the function of genes affected by a drought-rewatering cycle in the shoot apical meristem. Using microarray chips, differentially expressed genes (DEGs) and differentially methylated regions (DMRs) were identified for each water regime. The rewatering condition was associated with the highest variations of both gene expression and DNA methylation. Changes in methylation were observed particularly in the body of expressed genes and to a lesser extent in transposable elements. Together, DEGs and DMRs were significantly enriched in genes related to phytohormone metabolism or signaling pathways. Altogether, shoot apical meristem responses to changes in water availability involved coordinated variations in DNA methylation, as well as in gene expression, with a specific targeting of genes involved in hormone pathways, a factor that may enable phenotypic plasticity.

Non-cellulosic polysaccharide distribution during G-layer formation in poplar tension wood fibers: abundance of rhamnogalacturonan I and arabinogalactan proteins but no evidence of xyloglucan.

MAIN CONCLUSION: RG-I and AGP, but not XG, are associated to the building of the peculiar mechanical properties of tension wood. Hardwood trees produce tension wood (TW) with specific mechanical properties to cope with environmental cues. Poplar TW fibers have an additional cell wall layer, the G-layer responsible for TW mechanical properties. We investigated, in two poplar hybrid species, the molecules potentially involved in the building of TW mechanical properties. First, we evaluated the distribution of the different classes of non-cellulosic polysaccharides during xylem fiber differentiation, using immunolocalization. In parallel, G-layers were isolated and their polysaccharide composition determined. These complementary approaches provided information on the occurrence of non-cellulosic polysaccharides during G-fiber differentiation. We found no evidence of the presence of xyloglucan (XG) in poplar G-layers, whereas arabinogalactan proteins (AGP) and rhamnogalacturonan type I pectins (RG-I) were abundant, with an apparent progressive loss of RG-I side chains during G-layer maturation. Similarly, the intensity of immunolabeling signals specific for glucomannans and glucuronoxylans varies during G-layer maturation. RG-I and AGP are best candidate matrix components to be responsible for TW mechanical properties.

Expression analysis of LIM gene family in poplar, toward an updated phylogenetic classification.

BACKGROUND: Plant LIM domain proteins may act as transcriptional activators of lignin biosynthesis and/or as actin binding and bundling proteins. Plant LIM genes have evolved in phylogenetic subgroups differing in their expression profiles: in the whole plant or specifically in pollen. However, several poplar PtLIM genes belong to uncharacterized monophyletic subgroups and the expression patterns of the LIM gene family in a woody plant have not been studied. FINDINGS: In this work, the expression pattern of the twelve duplicated poplar PtLIM genes has been investigated by semi quantitative RT-PCR in different vegetative and reproductive tissues. As in other plant species, poplar PtLIM genes were widely expressed in the tree or in particular tissues. Especially, PtXLIM1a, PtXLIM1b and PtWLIM1b genes were preferentially expressed in the secondary xylem, suggesting a specific function in wood formation. Moreover, the expression of these genes and of the PtPLIM2a gene was increased in tension wood. Western-blot analysis confirmed the preferential expression of PtXLIM1a protein during xylem differentiation and tension wood formation. Genes classified within the pollen specific PLIM2 and PLIM2-like subgroups were all strongly expressed in pollen but also in cottony hairs. Interestingly, pairs of duplicated PtLIM genes exhibited different expression patterns indicating subfunctionalisations in specific tissues. CONCLUSIONS: The strong expression of several LIM genes in cottony hairs and germinating pollen, as well as in xylem fibers suggests an involvement of plant LIM domain proteins in the control of cell expansion. Comparisons of expression profiles of poplar LIM genes with the published functions of closely related plant LIM genes suggest conserved functions in the areas of lignin biosynthesis, pollen tube growth and mechanical stress response. Based on these results, we propose a novel nomenclature of poplar LIM domain proteins.

Microgenomic analysis reveals cell type-specific gene expression patterns between ray and fusiform initials within the cambial meristem of Populus.

The vascular cambium is the meristem in trees that produce wood. This meristem consists of two types of neighbouring initials: fusiform cambial cells (FCCs), which give rise to the axial cell system (i.e. fibres and vessel elements), and ray cambial cells (RCCs), which give rise to rays. There is little molecular information on the mechanisms whereby the differing characteristics of these neighbouring cells are maintained. A microgenomic approach was adopted in which the transcriptomes of FCCs and RCCs dissected out from the cambial meristem of poplar (Populus trichocarpa x Populus deltoïdes var. Boelare) were analysed, and a transcriptional database for these two cell types established. Photosynthesis genes were overrepresented in RCCs, providing molecular support for the presence of photosynthetic systems in rays. Genes that putatively encode transporters (vesicle, lipid and metal ion transporters and aquaporins) in RCCs were also identified. In addition, many cell wall-related genes showed cell type-specific expression patterns. Notably, genes involved in pectin metabolism and xyloglucan metabolism were overrepresented in RCCs and FCCs, respectively. The results demonstrate the use of microgenomics to reveal differences in biological processes in neighbouring meristematic cells, and to identify key genes involved in these processes.

Genome-wide analysis of LIM gene family in Populus trichocarpa, Arabidopsis thaliana, and Oryza sativa.

In Eukaryotes, LIM proteins act as developmental regulators in basic cellular processes such as regulating the transcription or organizing the cytoskeleton. The LIM domain protein family in plants has mainly been studied in sunflower and tobacco plants, where several of its members exhibit a specific pattern of expression in pollen. In this paper, we finely characterized in poplar six transcripts encoding these proteins. In Populus trichocarpa genome, the 12 LIM gene models identified all appear to be duplicated genes. In addition, we describe several new LIM domain proteins deduced from Arabidopsis and rice genomes, raising the number of LIM gene models to six for both species. Plant LIM genes have a core structure of four introns with highly conserved coding regions. We also identified new LIM domain proteins in several other species, and a phylogenetic analysis of plant LIM proteins reveals that they have undergone one or several duplication events during the evolution. We gathered several LIM protein members within new monophyletic groups. We propose to classify the plant LIM proteins into four groups: alphaLIM1, betaLIM1, gammaLIM2, and deltaLIM2, subdivided according to their specificity to a taxonomic class and/or to their tissue-specific expression. Our investigation of the structure of the LIM domain proteins revealed that they contain many conserved motifs potentially involved in their function.

Poplar genes encoding fasciclin-like arabinogalactan proteins are highly expressed in tension wood.

• Fifteen poplar cDNA encoding fasciclin-like arabinogalactan proteins (PopFLAs) were finely characterized, whereas the presence of arabinogalactan proteins (AGPs) was globally assessed during wood formation. • PopFLAs transcript accumulation was analysed through EST distribution in cDNA libraries, semi-quantitative RT-PCR, microarray experiment and Northern blot analysis. Similarly, AGPs contents were globally quantified by rocket electrophoresis. AGPs accumulation was further examined by Western blotting and immunocytolocalization. • Ten PopFLAs were specifically expressed in tension wood (TW) and not expressed in the cambial zone. Rocket electrophoresis revealed important AGPs accumulation in TW xylem. An anti-AGPs specific antibody recognized two proteins preferentially present in the cell wall-bound fraction from TW. Immunocytochemistry revealed a strong labelling close to the inner part of the G-layer of TW fibres. • PopFLAs are expressed in xylem and many are up-regulated in TW. It is suggested that some PopFLAs accumulating at the inner side of the G-layer may have a specific function in the building of this layer. PopFLAs expression may therefore be linked to the specific mechanical properties of TW.